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A practical parametrization for line shapes of near-threshold states 
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Numerous quarkonium(like) states lying near S- wave thresholds are observed experimentally.We 
propose a self-consistent approach to these near-threshold states compatible with unitarity and ana- 
lyticity.The underlying coupled-channel system includes a bare pole and an arbitrary number of elas¬ 
tic and inelastic channels treated fully nonperturbatively.The resulting analytical parametrization 
is ideally suited for a combined analysis of the data available in various channels that is exemplified 
by an excellent overall description of the data for the charged ^6(10610) and ^(10650) states. 
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INTRODUCTION 

At present there is no doubt that QCD is the true the¬ 
ory of strong interactions, at least at the energy scale 
presently accessible for experimental investigations. One 
of the remarkable features of QCD is the prediction of 
the existence of multiconstituent states, with a structure 
more complex than just quark-antiquark or three-quark 
configurations, which are conventionally referred to as 
“exotic” hadrons. Experimental searches and theoretical 
studies of such exotic states constitute an important tool 
in investigations of nature. Since the discovery of the 
charmonium(like jQ] state Y(3872) in 2003 0 , numerous 
experiments continue to deliver intriguing data on other 
charmonium(like) and bottomonium(like) states lying 
above the respective open-flavor thresholds. Although 
for most of these states it is not possible at present to 
make definite conclusions concerning their nature, some 
of these states share an important feature, namely they 
reside in the vicinity of strong S'-wave thresholds and 
they are seen in both open-flavor (elastic) and hidden- 
flavor (inelastic) final states. Paradigmatic examples of 
such states are the X(3872) near the D°D*° threshold, 
the Zb( 10610) and Zb{ 10650) near the B^B* thresh¬ 
olds, and the if c (3900) and Z c (4020)/Z c (4025) near the 
D^D* thresholds. Since vast and detailed information is 
becoming available from existing experiments, and even 
more precise data are expected from future high-statistics 
and high-precision experiments for states that are 
already known (see [6|, |7| for recent reviews) as well as 
for ones that are new and as yet unobserved, adequate 
theoretical tools for the data analysis are urgently called 


1 We refer to hadrons as to “charmonium(like)” or “bottomo- 

nium(like)” if they contain cc or bb quark-antiquark pair, re¬ 
spectively, however may have extra constituents like light-quark 
pairs. 


for. The aim of this Letter is to propose such a tool that 
is especially useful in describing near-threshold phenom¬ 
ena. 


The traditional way to perform an analysis of the ex¬ 
perimental data is by using of individual Breit-Wigner 
distributions for each peak combined with suitable back¬ 
ground functions. However, such an approach provides 
only limited information on the states studied, since the 
Breit-Wigner parameters are reaction-dependent and the 
naive algebraic sum of the Breit-Wigner distributions vi¬ 
olates unitary. In addition, by analyzing each reaction 
channel individually, one does not exploit the full in¬ 
formation content provided by the measurements. The 
approach proposed in this Letter provides an important 
link between various models and first-principles calcula¬ 
tions in QCD (for example, lattice simulations) from one 
side to the experimental data on the other side. To this 
end, we build a model-independent parametrization for 
near-threshold phenomena consistent with requirements 
from unitarity and analyticity. The formulas derived al¬ 
low one to perform a simultaneous analysis of the en¬ 
tire bulk of data for all decay channels for given near¬ 
threshold states(s). The resulting parametrization in¬ 
cludes in a fully nonperturbative way a bare pole and an 
arbitrary number of elastic and inelastic channels. With 
the help of not-very-restrictive and phenomenologically 
justifiable assumptions the formulas can be solved ana¬ 
lytically, which makes them as ideal for data analysis. 
The parameters of the final expressions are renormalized 
quantities with a direct physical meaning. The suggested 
parametrization is, therefore, expected to have a broad 
impact on the analysis of experimental data and to pro¬ 
vide important insights into the phenomenology of the 
strong interactions. 
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PARAMETRIZATION FOR THE LINE SHAPES 

We consider a coupled-channel approach based on the 
Lippmann-Schwinger equation (LSE) for the t matrix t, 


The omission of rescatterings within the inelastic chan¬ 
nels allows us to disentangle the latter from the elastic 
channels and from the pole term. We define the poten¬ 
tials 


t = V- VSt, 


(1) 


where S denotes the free propagator in the corresponding 
channel. The potential 


Pole 


0=1,JV e 




v = 


( 0 fp(jp') fi(k) \ Pole 

fa{p) V a p(p,p') V a i(p, k) 


= 1, N e 

\fj( k ') Vjp(k',p') Vji(k',k)J j = TJ^ 


( 2 ) 


contains all possible types of interaction between the bare 
pole (labeled as “0”—for example, its position is M 0 ), 
the set of N e elastic open-flavor channels ( Qq){qQ ) (here 
Q and q denote a heavy and a light quark, respectively) 
labeled by Greek letters, and a set of N in inelastic hidden- 
flavor channels ( QQ){qq ) referred to by Latin letters. 

In order to proceed with the analytic solution, we make 
a few simplifying assumptions. In general there are good 
reasons to neglect the direct interactions in the inelas¬ 
tic channels. For example, for the X(3872) transitions 
between the pJ/ip and u> J/ip channels are forbidden by 
the isospin conservation. In addition, since there are no 
light quarks inside the J/ip state, the direct potential 
for p(u})J/ip —> p(u})J/tp is also expected to be weak. 
Analogously, since there are no light quarks in the heavy 
quarkonia T (nS) and hb(mP), their interaction with pi- 
ons is expected to be weak, with obvious relevance for 
the states. Indeed, effective field theory estimates 
Q and lattice calculations |sj give very small values for 
the scattering lengths of the pion scattered off the cc and 
bb quarkonia. We therefore set Vji(k',k ) = 0. Next, 
we assume a separable form of the clastic transition ver¬ 
ted, v a i{p,k) = Xot{p)Via{k), where the additional as¬ 
sumption was made that Xia is independent of i. In¬ 
deed, the transition of the open-flavor channels to the 
hidden-flavor channels demands the exchange of a heavy 
meson and, therefore, it is necessarily of short range for 
all inelastic channels. Without loss of generality we set 
Xa(p = 0) = 1. In addition, in a relatively narrow region 
near the elastic threshold(s) it is sufficient to parametrize 
the transition form factors as 


fa(p) = f a ,Xa(P ) = 1 ,Via.{k) = g ia \ k\ li , /,; (fc) = A* | fc | ^ , 

where f a , gi a , A * are constants and U is the angular mo¬ 
mentum in the *th channel. The elastic potential, v a p, is 
approximated by a constant matrix. 


2 A microscopic model for this interaction can be found, for exam¬ 
ple, in Ha If ill. 


VaO = 


v oo = ~Y, = -E A * J - 

i i 

^ ^ y = — fa ^ ^ 

i i 

v 0( 3 = = u - 

i i 

~ = Va P ~ 5Z dia-Jigifi, 




where the thin solid lines, broad solid lines and dashed 
lines denote heavy-light mesons, heavy mesons, and light 
mesons, respectively, and the double line denotes the pole 
term. The inelastic loop integral is 

J i = j\q\ 2h Si(q)d 3 q -A 1 ^ m th in/il n (krf li+1 , (3) 

where the real part is omitted since it only renormalizes 
parameters of the interaction; p,\ n , k™, m th i n are the 
reduced mass, the relative momentum, and the threshold 
in the zth inelastic channel, respectively. To disentangle 
the pole term from the elastic channels we define 

= v a p — Gaf3 — V a0 G 0 Vop, Vao = ^4o(l + Go Too); 

where Go = l/(Mo — M + Vqo — *0) while the inelastic 
“bubble” operator is 


Gap — 'y ] ^ — y ' 9iaJi9if3 ■ (4) 

i i 

We arrive, therefore, at a pair of decoupled LSE 

t a p = yf P - E ^^p, t a o = v:$ - e 

r 7 p . (5) 

J a — I S a {^p)d p — (2tt)“^ q (k q T ika) = Pa T Hui 

with p, a and k a being the reduced mass and the rela¬ 
tive momentum in the ct’s elastic channel, respectively, 
k a = yj2p a (M — iTith a ) + where rrith a is the position 
of the ath elastic threshold. We reduced the entire prob¬ 
lem to Eqs. ©• Thus, independent of the number of 
inelastic channels the solution of these implies only the 
inversion of matrices as small as N e x N e , where typically 
N e = 2 (cf. the explicit example below). The transitions 
to inelastic channels follow from the solutions to Eqs. © 
straightforwardly, without the need to solve another scat¬ 
tering equation. Therefore, the proposed approach dras¬ 
tically simplifies the combined analysis of experimental 
data. In particular, adding a further inelastic channel 
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changes the final expressions only marginally. Since the 
approach is based on a LSE, unitarity is preserved auto¬ 
matically and all imaginary parts are linked to observable 
rates. 

In order to solve Eqs. © we proceed stepwise, analo¬ 
gous to the two-potential formalism Hl3|. Our starting 
point is a convenient parametrization for t v , the solution 
of the LSE t v = v — vSt v , where v is the direct inter¬ 
action potential in the elastic channels. The coupling to 
the inelastic channels is then switched on, and a LSE for 
the potential w = v — G, where the matrix G was defined 
in Eq. ©, is solved with the result 

r = t v + ^[g-G~ 1 }- 1 ^, (6) 

where the dressed vertices and the matrix Q are 

lfa.0 da0Xot ^0.0^0, ^fa.0 ~ da0Xa dofot0i 


Ga0 — Ja^al.3 — IpapJp — docpJa da.t 0 • 


( 7 ) 


Finally, when the coupling to the pole term is included 
as well, the formalism of [1J, ll5| can be used to yield 


tad — t 


af> M - M 0 + Go ’ 


taO — 


M — M 0 
M-Mo+S<> 


where 


4>a — baO ^ ' tgfiJpVp Q, 4>a — b()a ^ ' Vp0j0tff a , 

0 0 
Go — ^ Ji T y ( VftctJa&a = ^ ^ oYa. 0- 

i a i a 

The t matrix t a i is fully determined by t a p and t a o, 


taO A 
M - M 0 


-y^t af) R 0 g ip \(kr) li . (8) 


Since our knowledge of most resonance properties 
comes from production experiments, we build the pro¬ 
duction amplitude in the elastic or inelastic channel x 
as 


M x 





(9) 


where it was assumed that the production proceeds 
through the N e pointlike elastic sources T a . We also 
assumed that the elastic t matrix possesses poles near 
threshold(s) and, therefore, the Born term in the elas¬ 
tic amplitude was neglected. The differential production 
rate can be obtained by integrating the standard expres¬ 
sion for the three-body decay 16| in the invariant mass 
m| 3 , neglecting the FSI with particle 3. Then 


dBi x 

dM 


\M x \ 2 p 3 k x 

327r 3 M t 2 ot r to t ’ 


M = TO 12 = \fs, 


( 10 ) 


where k x is the c.m. momentum of particles 1 and 2. 
The allowed parameter range for M is given by M m i n = 
mi + m 2 and M max = M to t - m 3 . 

It is convenient to introduce new parameters A = JF\ 
and £ a = T a /T \, where the sources T a were redefined to 
absorb the slow function of energy R a and the constant 
factors from Eq. m- Since for all elastic channels the 
range of forces is described by the same physics, it is 
natural to use R a = ( 2ir) 2 p a n . Then the elastic and 
inelastic differential rates 


dBr e a 

dM 

dBr™ 

dM 


= A |y 

0 


2 

p 3 k a , 



PoK 


( 11 ) 

( 12 ) 


are described by the following set of parameters: 


A, Sa, fa, ^i, Qia, Mq, K, t V . 


(13) 


LINE SHAPES OF THE ^(10610) AND Z 6 (10650) 
STATES 


As an application for our approach, we consider l -1 
Zb( 10610) and Zb{ 10650) states residing near the BB* 
and B*B* thresholds, respectively, that are produced in 
T(5 S) decays T(55) —> nZb and that are seen in seven 
decay channels: Zb —> BB* , B*B* , nT(nS) (n = 1,2,3), 
and TrhbimP) (m = 1,2) [17l.ll8l|. The quantum numbers 
of the final quarkonia fix the angular momenta of the 
inelastic final states in Eq. (fl2l) to l = 0 for all nY(nS) 
final states and to l = 1 for the nhb(mP) final states. 

The fact that the 6-quark mass mb 3> Aqcd allows us 
to use heavy-quark spin symmetry (HQSS) to reduce the 
number of parameters. If the wave functions of negative- 
parity heavy-light B mesons are taken in the form (the 
charge conjugation is defined as CM = M) B = 0~g, 
B = 0^-, B* = l“g, and B* = —1^, then the Fierz 
transformation yields the following l -1 combinations: 


I BB*) 1+ - 
| B*B*) 1+ 


1 

1 

71 


(7b ® 0 m -) S= l + (0 b5 <g> lgg)s=l 

(K ® °7)s=i - (QfTb ® i««)s=i 


which imply that 0,0 


5[irT(nS)][B*B*] _ 1 9[-Kh b (mP)][B» B-] _ 1 

9[TrT{nS)][BB *] 9[-Kh b (mP)][BB *] 


(14) 


where the total angular momentum of the light-quark 
contribution in the latter case is to be provided by one 
unit of angular momentum that is explicitly accounted 
for in Eq. (TH?1) . Once the elastic channels BB* and B*B* 
are produced in the decays of the Y(5 S), the ratio of the 






















Arbitrary units 
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FIG. 1. (a)-(d) Fitted line shapes of the Z{,(10610) and ^(10650) states in the B^B* and Tvhb(mP) (m = 1,2) channels and 
(e) the predicted line shape in the 7 tY( 1S) channel [plots for 7rT(2S) and 7rY(3S) look similar and are not shown). Vertical 
lines indicate BB* and B*B* thresholds. Fit results are given in Eq. (IlSl- 


sources £ is subject to the same heavy-quark constraint, 

^ gfrT(5g)][B*B*] 1 

9[nT(5S)][BB*] 

In the same limit, the direct interaction in the B^B* 
system can be parametrized in terms of only two pa¬ 
rameters, 7 S and 7 1 , which are related to the contact 
potentials used in [ 21 | as 7 J 1 = (27r) 2 /r(Co a + Co b) and 
7t _1 = (27r) 2 ^(C 0a ^ 3C 0 h) in 1 « M 2 = n)- Then ^ 

+v = 1 1 ( hila + It) + ik 2 lilt-Is) \ 

(27 i) 2 n Det V lilt- Is) lils + It) +iki ) ’ 

with Det = 7 s 7 t - kik 2 + (*/ 2 )(Ys + lt)ih + k 2 ). 

The bare pole is included in the formalism in order to 
provide more flexibility in the fitting process- in partic¬ 
ular, it allows one to have two poles near the threshold 
even in the single-channel case. However, it should be 
omitted if its presence is not requested by the data. Thus, 
since we get a very good fit even without the bare pole, 
we refrain from its inclusion, thus setting f a = Aj = 0 
and Mq —> 00 in all formulas. As an experimental input 
we use 

• background-subtracted and efficiency-corrected 
distributions in M for the B^B* and Trhb(nP) 
channels [T3, [13 with floating normalization in 
each channel; 

• ratios of total branching fractions, Br^./Br^,^,, 

Br- n /Br^,^, {T 3 , [18, 0 , HU, where index i 

runs over all five inelastic channels ttT (nS) and 
7 r hb{mP). 

We do not use the information on the zj'- 1 line shapes 
in the 7 tT (nS) channels, since in the one-dimensional 
fit it is not possible to correctly take into account the 
interference with the nonresonant continuum, which is 
significant in the Y(5 S) —> n + TT~T(nS) transitions. In¬ 
clusion of the 7 tY (nS) line shapes would require a mul¬ 
tidimensional analysis, which is beyond the scope of this 
Letter; however, it is straightforward from the theoreti¬ 
cal point of view. In order to come to a converging fit 


we are, therefore, forced to impose that flWr(nS)][.B.B*] = 
- 9[nr(nS)][B*B *] for n = 1,2,3, as given by Eq. (TO). 
Meanwhile, we leave 3[^t(5S)][bb*] and 9[nT(5S)]{B*B*] 
as well as gfrh b (mP)][BB*] an d 9[*h b (mP)][B*B*] uncon¬ 
strained. The line shapes in the 7rY (nS) channels come 
out as our prediction. To take into account the experi¬ 
mental resolution, we convolve all the distributions with 
a Gaussian function with a = 6 MeV. Results of the 
simultaneous fit are shown in Figs. [U(a) [TJd). The de¬ 
veloped parametrization provides a very good description 
of the experimental data, with a confidence level of 76%. 
Predicted line shapes in the nT(nS) channels look rea¬ 
sonably similar to the experimental data [ 23 }; an example 
of such a distribution is shown in Fig. [1(e). It turns out 
that parameter k, defined in Eq. ©, is practically un¬ 
constrained by the fit; thus, we fix it to 1 GeV. From the 
fit we find 


£ = g[7rT(5S)1[5tjB ' >1 = —0.84 ± 0.05, 

5[7rT(5ST)][BB*] 

9{nh b (lP)]{B*B*] = 2 4± 0 6 

9[irh b (lP)][BB*] 

9[iTh b (2P)]{B*B «] = 2 4 ± 0 7 
9[nh b (2P)][BB*] 


(16) 


Deviations of these parameters from the predictions of 
HQSS—see Eqs. (fill) and (fl5l) -might be explained by 
the close proximity of the t matrix poles to the thresholds, 
which can result in an enhancement of the small explicit 
symmetry violation caused by Aqcd /mb ^ 0 01 • Are 
other source of the deviation of £ from the prediction of 
the HQSS may stem from a D -wave bb component as well 
as from possible non-66 components of the T(5S') wave 
function (cf. the discussion in [0). The importance of 
the HQSS-breaking contributions for the proper descrip¬ 
tion of the z)' 1 line shapes was also stressed in 261. It 
should be noticed that preliminary Belle data on the BB* 
and B*B* channels are used in the current analysis; fit 
results could change for the final experimental data. The 
inclusion of the information on the nT ( nS) line shapes in 
a future multidimensional analysis will help to improve 
the accuracy of the determination of the model param- 

























































eters and will allow for drawing more firm conclusions 
about the underlying physics of the spectacular near¬ 
threshold phenomena called Zb{ 10610) and i?b(10650). 
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CONCLUSIONS 


In this Letter we proposed a practical parametriza- 
tion for the line shapes of the near-threshold state(s). 
Since the approach is based on the LSE for the coupled- 
channel problem, unitarity and analyticity constraints for 
the t matrix are fulfilled automatically. This guarantees 
that all imaginary parts are included in a self-consistent 
way. Since there are good reasons to neglect direct in¬ 
teractions within the inelastic channels, at least for the 
systems discussed here, the inelastic channels enter the 
expressions additively which; this makes it particularly 
easy to extend the inelastic basis. While additional ef¬ 
fects such as finite widths of the constituents and the 
FSI with the spectator may also play a role and should 
be included on top of the interactions considered in this 
Letter; nevertheless, the gross features of the coupled- 
channel problem are captured by the presented model 
and the parametrization based on it is expected to be 
realistic. Finally, we demonstrate the power of the sug¬ 
gested parametrization by the fit to the line shapes for 
the bottomoniumlike states Zb{ 10610) and Zb{ 10650), for 
which we obtain a very good description. 
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